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Abstract 

We review the tidal capture mechanism for binary formation, an important process 
in globular cluster cores and perhaps open cluster cores. Tidal capture binaries 
may be the precursors for some of the low-mass X-ray binaries observed in abun- 
dance in globular clusters. They may also play an important role in globular cluster 
dynamics. We summarize the chaos model for tidal interaction (Mardling 1995a), 
and discuss how this affects our understanding of the circularization process which 
follows capture. 

1.1 Introduction 

Various theories for binary star formation have been suggested (some more suc- 
cessful than others), including fission (Durisen & Tohline 1985), fragmentation 
(Boss 1985), conucleation (Shu et al. 1987) and dynamical formation and tidal 
capture (Clarke, this volume), all of which are supposed to occur around the time 
of star formation. Dynamical formation (see next section) depends only on the 
masses of the objects and requires high protostellar densities to be effective, as are 
found in protostellar clusters. This is also true of stellar tidal capture, a process 
which becomes viable in dense stellar environments such as globular cluster cores 
(and perhaps open cluster cores). We will describe tidal capture in detail in the 
next section, but suffice to say here that it involves transferring the excess energy 
of unbound orbital motion to the tides of two stars which pass each other at a 
distance of a few stellar radii, so that a bound system (a binary) can result. 

1 Permanent address: Mathematics Department, Monash University, Melbourne, Australia, 
r.mardling@sci.monash.edu . au . 
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The concept of tidal capture occurred to Fabian, Pringle & Rees in 1975 af- 
ter several X-ray sources were discovered which appeared to be associated with 
globular clusters. It was inspired by the fact that the X-ray luminosity-to-mass 
ratio was more than 100 times greater in globular clusters than for the galaxy as a 
whole (Katz 1975), so that the mechanism responsible for the X-ray sources must 
be peculiar to globular clusters (but see Mardling 1996). They showed that it was 
possible for capture to occur in a globular cluster (as opposed to only collision or 
flyby) and proposed that the X-ray sources were accreting neutron stars. Press 
& Teukolsky (1977) followed up this work by calculating the energy transfer to 
the tides of a polytrope of index 3 during periastron passage of a parabolic orbit 
(see also Lee & Ostriker 1986 and Giersz 1986). At around the same time, it was 
suggested by Gunn & Griffin (1979) that globular clusters were devoid of binaries. 
Since it was becoming clear that binaries could provide a source of energy for core 
support against collapse, 2 tidal capture binaries took on this role as well as the 
role of progenitors for the X-ray sources. 

Later McMillan, McDermott & Taam (1987) suggested that the orbital circu- 
larization process following capture destroys the binary (see also Ray, Kembhavi 
& Antia 1987). They argued that energy transfer at periastron is one way (from 
orbit to tides) so that the binary must circularize very quickly (in as little as 10 
yr). They estimated a dissipation timescale of 10 4 — 10 6 yr so that by the time 
of circularization, a system has dissipated only a small fraction of the total tidal 
energy. The tidal energy will be equal to the binding energy of the binary since 
the binary orbit at capture is approximately parabolic. Assuming that the orbital 
angular momentum remains constant ( ie. assuming there is no angular momentum 
loss via mass loss, gravitational radiation or magnetic breaking), the final orbital 
separation will be twice the periastron separation at capture. For an equal-mass 
system with a capture periastron separation of 3 stellar radii, the tidal energy will 
be 20% of the internal binding energy of the polytrope. The stars will respond 
by expanding; since they are so close initially, a common envelope phase will be 
entered into and the binary may be destroyed. 

But this didn't leave globular clusters devoid of binaries; finally some giants 
were observed in the outer regions of a few clusters, from which a binary fraction 
of about 10% was inferred (Pryor et al. 1989). 

In this way tidal capture fell from favour as an important process in globular 
clusters. But it still remained to explain the origin of the X-ray sources, long 
agreed to be low-mass X-ray binaries (LMXBs - neutron stars accreting from 
low mass main sequence or giant companions). Hut, Murphy & Verbunt (1991) 
suggested that such binaries were the result of clean exchanges of neutron stars 
into primordial binaries. However, this process tends to leave binaries with periods 
considerably longer than those observed (Mardling 1995b). Since the mass ratio 
of a neutron star to a globular cluster main sequence star (the turnoff mass is 
~ 0.8M Q ) is such that a common envelope phase will not be entered into, the 
only obvious way to shrink the orbit within a Hubble time is through encounters 

2 Indeed it was shown using Fokker-Planck simulations that the tidal capture binaries produced 
by a cluster are capable of halting core collapse and could even fuel reexpansion (Statler, Ostriker 
& Cohn 1986). 
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Figure 1.1: Dynamical formation of a binary. Here Etot and Jtot are the 
total energy and angular momentum respectively, Eij is the orbital binding 
energy of the pair and Jij is the angular momentum about the centre of 

mass of the pair. 



with other stars. It remains to be seen whether this process is efficient enough to 
produce the LMXBs observed (Davies 1995). 

And it still remained to understand what actually happens following capture. 
We will address this question in the following section in which we describe in more 
detail the tidal capture process. Section 2 summarizes the chaos model for tidal 
interaction (Mardling 1995a), while section 3 applies this model to capture orbits. 
Section 4 discusses dissipation and long-term evolution, while the final section 
presents a discussion. 

1.2 Tidal Capture 

Consider an isolated system of two point masses moving in a hyperbolic orbit 
relative to each other. It is not possible for such a system to produce a binary; 
the relative orbit is fixed along with the the total energy and angular momentum. 
Now introduce a third body into the system. Fig. 1.1 illustrates how it is possible 
to produce a binary from the original pair by transferring some of its orbital 
energy and angular momentum to the third body. This process of forming binaries 
dynamically occurs in dense clusters; in fact it is inevitable in a cluster consisting 
only of point masses and which contains no binaries initially (Aarseth 1974). 

Now return to two bodies but allow one of them to be finite in size and non- 
rigid. Let their relative speed at infinite separation be so that the total energy 
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Figure 1.2: Tidal capture formation of a binary. .Eorb and E osc are the 
orbital binding energy and the tidal or oscillation energy respectively, and e is 
the orbital eccentricity. 



of the system is ^/w^,, where \jl = M\Mij{M\ + M^) is the reduced mass with M\ 
and Mi the masses of the particles. Now it is possible for the non-rigid body to 
absorb some of the energy of the system; if it can manage to absorb at least \^v'l , 
the orbital binding energy of the pair will become negative and a binary will have 
been formed. This is the process of tidal capture (see fig. 1.2). Now let us examine 
the energy transfer process. It is easier to imagine the extended object as being 
composed of N particles rather than a continuous fluid. Each particle experiences 
a force from the companion object, and this force will vary from particle to particle 
depending on its distance from the companion. At the same time, each particle 
experiences the gravitational force due to all the other particles, as well as a 
pressure force. As the companion draws near, the particles are pulled towards it 
and the fluid body assumes a distorted shape. But the restoring forces of pressure 
and self-gravity respond and an oscillation is set up. This process is commonly 
referred to as dissipative because in all but ideal circumstances, the oscillations 
will damp via viscous processes within the star. We refer to the oscillation that is 
set up as the tides of the star. 

Just as in the case of a violin string or a drum, the oscillation can be decom- 
posed into the object's normal modes of vibration. A string being 1-dimensional 
has one set of eigenmodes which are the circular functions. A drum has two sets; 
the radial eigenfunctions are Bessel functions, while the azimuthal eigenfunctions 
are again the trigonometric functions. The normal modes of a star will be made 
up of three sets of eigenfunctions. Given an equation of state for the fluid the 
radial eigenfunctions can be determined. These will come from a system of equa- 
tions which resemble a Sturm-Liouville system. The polar angle dependence is via 
Legendre functions and the azimuthal dependence is trigonometric, these two to- 
gether forming the spherical harmonics (Jackson 1975). The modes exited in a star 
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Figure 1.3: Maximum periastron separation for capture, i?™ ax (in units of 
the stellar radius), versus mass ratio q for point mass/polytrope models with 
indices n = 1.5 and 3. Here q is the mass ratio of the point mass to the 
polytrope and «oo is the relative velocity at infinity measured in km s _1 . 



by a companion are called the non-radial modes of vibration because they have 
angular dependence. Purely radial modes require a radially symmetric excitation 
mechanism as in the case of the Cepheid variables. 

Whether or not a capture results in a stable binary (at least in the short term) 
will depend on how close the two stars come as well as the system's total energy, 
all of which must be absorbed by the tides. If the total energy of the system 
is excessive, the stars will essentially need to collide in order to transfer enough 
energy to create a bound system. If the stars avoid collision and manage to form 
a binary, they may still tidally disrupt on a subsequent periastron passage. 



1.2.1 Capture Cross-Section 

How close do stars need to come for a capture to occur? For velocity dispersions 
typical of open and globular clusters (1 km s _1 and 10 km s _1 respectively), their 
separation at closest approach must be at most a few stellar radii (fig. 1.3). But the 
process is significantly enhanced by the gravitational focussing of the two bodies. 
Fig. 1.4 illustrates how this works with two point masses. Let Rq be the perpen- 
dicular distance between the two bodies at infinity (the impact parameter) and 
R p (p for periastron) be the separation at closest approach, and let v p = R p (f p be 
the velocity of the companion relative to the first star, with tp p being the angular 
velocity at periastron. Conservation of angular momentum gives 

.Ro^oo = RpVp (1.1) 
while the total energy of the system is given by 

1 2 1 2 GM 1 M 2 

1 HV 00 = 1 nv p . (1.2) 
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Table 1.1: The effectiveness of gravitational focussing 



Voo (km s -1 ) 


-ftp 


Ro 


10 


3.2r* 


1 AU 


1 


4.1r, 


12.6 AU ~ 1.3i? Sa turn 



These two equations together give an expression for the capture cross-section, a: 
a = nR 2 = ttR 2 p (l + (v esc / Voo ) 2 } , (1.3) 

where v esc = (2G(M 1 + M 2 )/R p ) 1 / 2 is the "escape velocity" of one star relative 
to the other at a distance R p . Since this is generally much greater than Voo for 
globular and open clusters, 3 the term involving it dominates and we have the 
result that a ~ R p . In contrast, the velocity dispersions in galactic nuclei are 
much higher so that the capture cross-section is not dominated by gravitational 
focussing in this case. 

Table 1.1 illustrates the effectiveness of gravitational focussing for the capture 
of a neutron star (mass 1.4M©) by a main sequence star of mass O.6M and 
radius 0.6i? Q . Here, i?™ ax is the maximum distance of closest approach for which 
a capture can occur, a quantity which depends on the mode structure of the star 
(in this case for a polytrope of index 1.5). This can be calculated using the model 
to be described in this paper, or using the analysis devised by Press and Teukolsky 
(1976) to calculate the energy transfer during one periastron passage of a parabolic 
orbit. 

1.2.2 Capture Rates 

How often does tidal capture occur between two species of stars? This will depend 
on the stellar densities of both species, their relative velocities, and the capture 
cross-section (which itself depends on the relative velocities). Imagine we have 
species 1 and 2, say, a population of neutron stars of mass Mi, and a population 

3 For example, the escape velocity at the solar surface is 619 km s -1 . 
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of main sequence stars of mass M 2 . For simplicity we can imagine that each 
neutron star-main sequence star pair is moving relative to each other with the 
same speed, in this case, the cluster velocity dispersion. If the stellar density of 
species 1 is N\ , then the number of captures of these stars by a star of species 2 
in time St will be ^(v^S^a, that is, the number of stars inside a volume element 
of side length v^St and base area er. If the stellar density of species 2 is N 2 , then 
the total number of captures that take place in a unit volume and unit time, that 
is the capture rate between species 1 and 2, is given by 

r 12 = Ni^v^a 

= 2ttGN 1 N 2 R p (M 1 + M 2 )/ Voo 

co N, N 2 Mi + M 2 Rf n 10 km s" 1 _ 3 
= 6 - 8 10>c-3l 04 pc-3 Mq -R^^— GyT PC • (L4) 

Using data from Davies & Benz (1995) for the globular cluster 47 Tuc, the density 
of neutron stars in the core is 900 pc~ 3 , the density of turnoff mass main sequence 
stars is 3 x 10 4 pc~ 3 and the velocity dispersion is 10 km s _1 so that the capture 
rate is 12 captures Gyr -1 pc~ 3 . 

1.2.3 Evolution Following Capture 

What really happens following capture? The energy transfer during the first peri- 
astron passage is well established, at least for non-rotating non-pulsating stars in 
parabolic orbits. But the energy transfer during subsequent passages will depend 
on whether or not the stars are still oscillating, which in turn will depend on the 
dissipation timescales of the stars compared with the orbital periods. For captures 
occurring with velocities at infinity typical of globular cluster (3-D) velocity dis- 
persions (around 10 km s _1 ), the orbital periods are less than the estimated linear 
dissipation timescale (see fig. 1.5), while for lower velocities (such as would be 
found in open clusters), the situation can reverse. On the other hand, for very en- 
ergetic tides nonlinear effects can become important and the dissipation timescale 
may be much shorter. But if the stars are still oscillating during a periastron 
passage, the direction of energy flow will depend on the phase of oscillation; it is 
possible for energy to be returned to the orbit (Mardling 1991, Kochanek 1992). 
Thus it is also possible for a newly formed binary to self-ionize if the system has 
not yet dissipated all the energy of the initial unbound orbit. But it can be shown 
that the majority of systems avoid this fate (Mardling 1995b). 

The energy transfer at periastron will also depend on the tidal energy already 
present. In fact as we will show in section 1.3, this energy transfer process is 
chaotic; the evolution of the orbit following capture depends sensitively on con- 
ditions at capture, at least until the system has dissipated an amount of energy 
which depends again on conditions at capture. 

1.2.4 Dynamical Calculations 

In order to study the dynamical evolution following capture, one must follow in 
detail the individual oscillations which occur on the dynamical timescale of the 
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Figure 1.5: The first orbital period following capture as a function of pe- 
riastron separation at capture for various velocities at infinity (measured in 
km s _1 ). The asymptotes indicate the maximum periastron separation for 
which capture is possible. 



star. But the orbital periods immediately after capture are extremely long. For 
example, for an equal-mass capture at 3 stellar radii with a velocity at infinity of 
10 km s _1 (corresponding to an eccentricity of 1.0008 before capture and 0.999 fol- 
lowing capture), the orbital period is 1.7x 10 5 oscillations of the lowest order (most 
energetic) mode. This fact has previously prevented the dynamical calculation of 
all but a few orbits following capture. 

Several numerical methods have been applied. Smoothed particle hydrody- 
namics (SPH) was used to study collisions and extremely close encounters (see 
for example Gingold & Monaghan 1980, Rasio & Shapiro 1991 and Benz & Hills 
1992). 

A method devised by Carter & Luminet (1985) which employs a polytropic 
version of Chandrasekhar's homogeneous ellipsoid was used by Kochanek (1992) 
to model the first few orbits following capture. 

A third method is the subject of the rest of this review. It is a modification 
of a normal mode analysis devised by Gingold & Monaghan (1980). The success 
of this method lies with the fact that the equations of motion are derived from a 
Lagrangian so that the total energy (and angular momentum) are conserved and 
energy is free to flow between the orbit and the tides. Details may be found in 
Mardling (1995a, b); here we present a summary. 

1.3 The Chaos Model for Tidal Interaction 

We now examine more formally the energy transfer process of the tidal interaction 
between two stars (here modelled as a point mass and a non-rotating polytrope) . It 
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is possible to construct a Lagrangian consisting of terms corresponding to a Kepler 
orbit, terms corresponding to a self-gravitating fluid and a term corresponding to 
the interaction between the two. Thus energy is free to flow either way between the 
orbit and the tides, and this interaction can be followed in detail. This calculation 
can be found in Gingold & Monaghan (1980) and Mardling (1995a); here we 
present an outline of how the method works. 4 

Given a polytropic equation of state, one then assumes small perturbations to 
the density and velocity field of the polytrope. The perturbation variables satisfy 
normal mode equations as discussed above. The Lagrangian is then expanded to 
second order in the perturbation variables and orthogonality conditions associated 
with the normal modes are used to perform any integrals over the polytrope. In 
this way, a Lagrangian is derived which depends on the orbital variables as well as 
the time-dependent mode amplitudes. One can then derive equations of motion 
for these variables and a conserved energy (and angular momentum) . Higher order 
perturbation expansions can also be done in this way 5 (Mardling 1991), and since 
the equations of motion are derived from a Lagrangian, one is always assured of 
the existence of a conserved energy. 

The equations of motion are in the following form: 

b k +Lot l b k = qCl 1) — [W , (1.5) 

= B^ K + q 2- ^ 6k 5R [r^ ) ' (L6) 

where 6k is the (complex) amplitude of the mode with mode numbers k = klm, 

oo oo / 

= ^2 ^2 ' ^ s frequency of mode k, tp is the true anomoly mea- 

k k=l 1=2 m=-l 

sured from the initial line of apses, R = |R| is the binary separation, q is the mass 
ratio of the point mass to the polytrope, and C k ^ and C k 2 ' are constants which 
depend on the internal structure of the polytrope. Eq. (1.5) represents a forced 
harmonic oscillator, the forcing being provided by the orbit, and eq. (1.6) is a 
Kepler orbit with a mode-orbit coupling term as a perturbation. The strength of 
the interaction between the orbit and the tides increases as the central condensa- 
tion of the polytrope decreases, or as the mass ratio is increased. The solution is 
constrained by an energy integral as well as an angular momentum integral: 

e = y& - - g £ c^K^ + (Mi + ^K), (i.7) 

k k 



4 Another way to set up this problem (Mardling 1991) is to regard the fluid as being composed 
of N particles and to consider the gravitational force between each of the particles, including the 
point mass, and the pressure force at each particle. The continuous limit is then taken. This way 
of approaching a self-gravitating problem can have many advantages, especially if the problem 
has odd geometry; in fact the method of SPH is set up in this way (without the continuous limit). 

6 Although one must be careful to avoid divergent series by allowing the frequencies of vibration 
to vary in time. 
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J = I^R 2 (p + J2 imC L 3)b ^l, (1.8) 

k 

where = /C^ and i = >/^T. The total energy B consists of the orbital 
binding energy (the first two terms) , the oscillation energy of the tides (the last two 
terms), and the interaction energy between the two. Note that angular momentum 
can be transferred to the tides in the absence of friction. It is commonly believed 
that angular momentum transfer is only possible when dissipation causes a tidal 
lag, that is, the tidal bulge of the equilibrium tide does not point in the direction of 
the companion. This in turn places a torque on the star, ie., angular momentum is 
transferred from the orbit to the tides. But the asymmetric nature of the dynamical 
tide allows angular momentum to be transferred in the absence of dissipation. 

The self-consistent system of equations (1.5) and (1.6) allows energy and an- 
gular momentum to be freely exchanged between the orbit and the tides, the 
consequences of which can be severe. 

1.3.1 Chaotic and Periodic Solutions 

Two types of solutions exist to the system of eqs (1.5) and (1.6). A convenient way 
of displaying these solutions is to plot the orbital eccentricity against the number 
of periastron passages. This is done in fig. 1.6 for three sets of initial conditions, 
with each orbit starting with an eccentricity of 0.8 (the polytropic index is 1.5). 
The widest orbit (a) has an initial periastron separation 6 of 3.2r*, where r* is the 
radius of the polytrope. The change in eccentricity is small from orbit to orbit, 
and a beating is evident in the enlargement shown at the bottom of the figure. 
Orbits (b) and (c) differ in their initial periastron separation, R p , by less than one 
part in 10 5 (R p = 2.9r* and R p = 2.90001r*). They exhibit extreme sensitivity to 
initial conditions indicating chaos. 

Why is the behaviour of chaotic and periodic systems so different? There 
appears to be no lower bound for the eccentricity of chaotic orbits (in fact there 
does exist such a lower bound; see section 1.4), while for periodic orbits, the 
maximum change in eccentricity is of the order of the change during one periastron 
passage. 7 Large changes in eccentricity imply very energetic tides which in turn 
imply a need to consider nonlinear mode couplings, something this model neglects 
(but see section 1.4). In contrast, the tidal interaction in a periodic orbit will 
generally raise small tides. The difference between the two behaviours can be 
understood mathematically from the point of view of stability. The author has 
devised a 3-D mapping which mimics the behaviour seen here. The system appears 
to follow the "intermittency route to chaos" (Pomeau & Manneville 1980). If one 
calculates the eigenvalues of the mapping, stability requires their absolute values 
to be unity (the system is non-dissipative). When a system goes unstable, ie. 
becomes chaotic, the eigenvalues wander arbitrarily far from unity. The details of 
this work will be published elsewhere. 

6 The orbit is actually started at apastron and the periastron separation is supplied which 
would be reached if the system consisted of two point masses. 

7 Except for resonant orbits - those for which the orbital period is nearly a multiple of the 
period of oscillation of the most energetic mode. 
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Figure 1.6: Chaos vs. periodicity. Curve (a) (shown magnified at the bottom 
of the figure) is typical of a periodic orbit for which the eccentricity varies 
little so that the tidal energy is small. In contrast, curves (b) and (c) display 
the extreme sensitivity to initial conditions typical of chaotic systems, with 
the eccentricity varying over a wide range and the tidal energy becoming very 
large. 



1.3.2 The Chaos Boundary 

The behaviour demonstrated in fig. 1.6 is typical of chaotic and periodic systems. 
We can discover the range of initial periastron separations and eccentricities for 
which chaotic motion exists by exploiting the extreme sensitivity to initial condi- 
tions. Fig. 1.7 plots the boundary between chaotic and periodic motion for equal- 
mass systems containing a polytrope of index 1.5. This boundary was derived by 
systematically running though (R p ,e) parameter space and comparing 50 perias- 
tron passages of two orbits whose initial eccentricity differs by 10~ 5 and whose 
initial tidal energy is zero. Naturally there exists the possibility that a system 
will appear periodic for more than 50 orbits and suddenly become chaotic. This is 
characteristic of systems which follow the intermittency route to chaos. Nonethe- 
less, this kind of behaviour is restricted to regions close to the chaos boundary so 
that although the boundary as plotted is not well defined, it serves the purpose of 
delineating the two types of behaviour. 

The position of the boundary depends on the mass ratio as well as the stellar 
structure. If we define the mass ratio to be the mass of the compact object to 
the mass of the polytrope, then increasing this ratio causes the chaos boundary 
to move up in (R p ,e) space. This reflects the fact that for a fixed separation, 




Figure 1.7: The chaos boundary. The polytropic index is 1.5 and the mass 
ratio of compact to extended star is unity. The boundary moves up if the mass 
ratio is increased, and down if the polytropic index (central condensation) is 
increased. The boundary for tidal disruption is underestimated by this model. 



the tidal force increases as the mass ratio is increased. This can also be seen in 
eqs. (1.5) and (1.6); the strength of the interaction is proportional to the mass 
ratio. Similarly, the chaos boundary moves down if the central condensation of 
the polytrope increases (equivalently, if the polytropic index increases) . For a fixed 
stellar radius and stellar separation, the energy the tides are capable of storing 
decreases with increasing central condensation. This reflects the mode structure 
which in turn reflects the mass distribution. 

The chaos boundary is not very sensitive to the following: 

1. The number of modes in the calculation. Including more than the I = 2 
/-mode (k = 1) has very little effect on the position of the chaos boundary. 

2. A finite radius companion. 

3. Including mode-mode interactions (ie. including higher order terms in the 
perturbation expansion) . 

Each case can be understood in terms of stability as discussed above; the eigen- 
values of the system are not much affected by these factors. 

An interesting feature in fig. 1.7 is the asymptote at e = 1: all (nondissipative) 
capture binaries are chaotic. The asymptotic nature of the chaos boundary can also 
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periastron passage 



Figure 1.8: A capture orbit with an initial eccentricity of unity and a perias- 
tron separation at capture of 3r* (no dissipation). Note the periodic feature 
around the 2000th orbit. At this point, the height £ of the tides is of the order 
of the stellar radius. 



be understood in terms of the mapping representation of this problem (Mardling, 
in preparation). As well there exists a second branch to the chaos boundary which 
corresponds to hyperbolic capture orbits. 

1.4 Capture Orbits 

We now apply the chaos model for tidal interaction to the problem of tidal capture. 
We start by considering a dissipationless model, a model which is highly unrealistic 
because dissipation plays a major role in the dynamical evolution following capture. 
Nonetheless, this enables us to examine certain properties of such systems. 

But first we must overcome the problem of computing very long period orbits 
while at the same time following the motion of the tides. We do this by recogniz- 
ing that the tides and the orbit essentially only interact for a fraction of a very 
long period orbit - that is for that section of the orbit near periastron. We thus 
calculate this section by numerically solving eqs (1.5) and (1.6) sufficiently far past 
periastron so that the total energy and angular momentum (eqs. (1.7) and (1.8)) 
are conserved to within some tolerance 8 after several thousand orbits. The rest of 
the orbit is solved for analytically (see Mardling 1995b). 

Figure 1.8 shows 3000 orbits of an equal-mass system in an initially parabolic 
orbit with an initial periastron separation of 3r* (the calculation includes up to 



8 To within 1 part in 10 4 in the example which follows. This corresponds to halting the 
numerical calculation at a stellar separation of about 40 stellar radii. 
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Figure 1.9: The evolution curve. Comparison with the curve along which the 
orbital angular momentum is constant indicates that a significant amount of 
angular momentum is transferred to the polytrope. Note the endpoint of the 
evolution curve which represents a minimum orbital eccentricity the system 
can achieve in the absence of friction. The finite width of this curve is mainly 
due to the influence of the 1 = 3 mode and relects its independence from the 
1 = 2 mode. 



the I = 4 mode for k = 1 only). The eccentricity varies over a wide range but 
there exists a lower limit, in this case e = 0.5. The quasiperiodic behaviour seen 
around the 2000th orbit indicates that the system has met a chaos boundary. 9 In 
the absence of dissipation, the eccentricity cannot drop below this value. 

Since the tides are extremely energetic (this calculation implies a tidal height of 
the order of the stellar radius when e = 0.5), a linear analysis is clearly inadequate. 
We argue though, that neglecting nonlinear mode-mode interactions is equivalent 
to neglecting a strong source of dissipation, at least for systems sufficiently far 
from tidal disruption. Since the orbit can only drive the low order modes (the tidal 
energy in a mode is proportional to R~( 2l+2 ^) and since the system is chaotic, it is 
impossible to distinguish between the present model and a model which includes 
mode-mode interactions (Mardling, in preparation). 

It is also instructive to plot R p against e and compare this evolution curve 
with the curve along which the orbital angular momentum is constant, the latter 
being given by R p (l + e) = J?°(l + eo), with (R°, eo) the initial orbital variables at 
capture. Fig. (1.9) shows this comparison along with an enlargement of a portion 
of the figure, showing that the evolution "curve" is actually of finite width and 
consists of points which wander up and down in a chaotic manner. The difference 
in the two curves reflects the angular momentum transferred to the tides, despite 
the system being frictionless. In fact, one can show that the angular momentum 
in a particular mode is proportional to the energy in that mode. 

The left-hand end of the evolution curve defines the minimum possible eccen- 



9 A non-zero tidal energy chaos boundary; see Mardling (1995a) for details. 
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Figure 1.10: The effect of dissipation: the evolution curve shrinks towards 
the point P after which the system becomes permanently periodic and evolves 
along the curve of constant angular momentum. 



tricity the system is capable of attaining (in the absence of dissipation), 10 which 
in turn limits the maximum possible tidal energy. In this example, the maximum 
possible tidal energy is 7% of the stellar internal binding energy, 11 quite a large 
fraction (although nonlinear dissipation is likely to severely limit this maximum 
- see the next section), while for wider systems, this maximum can be arbitrarily 
small. We will see in the next section that this has important implications for the 
survival of binaries after capture. 

1.5 Dissipation and the Long-Term Evolution 

We have neglected dissipation so far in order to demonstrate the chaotic nature 
of capture systems. Assuming there is no mass loss or transfer, and that the radii 
of the stars remains constant, it is possible to examine qualitatively the effect of 
dissipation on the orbital dynamics without knowing the details of the dissipa- 
tion mechanisms. Consider the effect on the evolution curve as the system loses 
energy. The maximum possible eccentricity will decrease as the system becomes 
more bound, and since this corresponds to zero tidal energy, the endpoint of the 
evolution curve must sit on the curve of constant orbital angular momentum, 
Jmh = const (fig. 1.10). At the same time, the minimum possible eccentricity will 
increase, a process which can be partially understood as follows (but see Mardling 

10 It is possible to calculate this minimum; see Mardling (1995b). 

11 This should be compared with a figure of 20% for complete circularization had the tidal 
energy not been restricted by the chaos boundary 
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Figure 1.11: A capture orbit with dissipation. 



1995b). Imagine that a system has evolved to the state in which the right hand 
end of the evolution curve as described above lies at the point where the chaos 
boundary Co intersects J or b = const (the point P in Figure 1.10). By definition, 
the right hand end point of the evolution curve corresponds to zero tidal energy. 
Since any system which starts to the left of the chaos boundary with zero tidal en- 
ergy is periodic in nature, it follows that the point P marks the end of the chaotic 
phase and that before this, the minimum eccentricity that the system can achieve 
at any time must increase. After the system has reached the point P and is no 
longer chaotic, the tides will be small and hence the dissipation rate will be greatly 
reduced. The system will then move slowly along J or b = const until the binary is 
circularized with a separation of twice the periastron separation at capture. 

This process can be modelled by introducing an artificial damping term into 
the equation of motion for the mode amplitudes. Eq. (1.5) is replaced by 

b k + 2C k b k + u>l l b k = qCl 1) — [Tr , (1.9) 

where oc 1/tu, with Tk being the damping timescale for mode k. In fact, 
when the tides are large, nonlinear damping due to mode-mode interactions will 
dominate linear damping, the latter being due to normal viscous processes. The 
damping timescale for mode k will thus be given by 1/rk = 1/T k + 1/t£ 1 , where r k 
and t£ 1 are the linear and nonlinear damping timescales respectively. Figures 1.11 
and 1.12 illustrate the circularization process again for an equal-mass system 
with initial periastron separation R p = 3r* and initial eccentricity e = 1, but 
with a nonlinear damping timescale of about 30 yr and a linear damping timescale 
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Figure 1.12: Chaotic and periodic circularization for the system shown in 
Figure 5. The times shown in days are the orbital periods at these points. 



of 1000 yr. The nonlinear damping timescale may well be shorter than this, but 
whenever it is longer than the orbital period, chaotic behaviour will persist. In this 
case, the chaotic phase lasts for about 200 years after which time the binary takes 
more than 5.5 million years to finish circularizing. This represents on the order 
of 2 x 10 9 orbital periods which would take around 1000 CPU years to run on a 
fast workstation. Instead this calculation was done using a method devised by the 
author for use in TV-body cluster calculations (see Aarseth 1996 for a description 
of the astrophysical processes included in this model). This method assumes an 
impulse approximation for energy transfer at periastron during the chaotic phase. 
The energy transfer depends on the tidal energy present as well as the oscillatory 
phase and uses a modified Press and Teukolsky (1977) analysis (one for which 
the similarity variable depends also on the orbital eccentricity). The periodic 
phase is calculated semi-analytically. A thorough description of the method will 
be published elsewhere. 

1.6 Discussion 

It is not clear how the stars will respond to the violent tides during the chaotic 
phase. The stars are so close during the circularization process that any expan- 
sion in response to tidal heating will completely change the evolution. Thus it is 
important to know where the tidal energy is deposited. Podsiadlowski (1996) has 
shown that the stellar radius is not much affected if the energy is deposited near 
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the core of the star. If the energy is deposited in the extreme outer layers, the star 
may lose part of these layers via a stellar wind resulting in expansion of the orbit, 
and this in turn will reduce the tidal interaction (although it should be noted that 
the chaos boundary will also be affected). On the other hand, if the energy is 
deposited in the intermediate layers, the effect could be devastating. Clearly the 
fate of tidal capture binaries depends on what happens during this phase, but one 
might expect a reasonable fraction (maybe 25%?) of binaries which survive the 
first few orbits following capture to go on and enter the quiescent phase intact. 

Of those binaries which do go on to circularize "normally" , how many binaries 
containing compact objects eventually become stable mass transfer systems? The 
mass ratio of a neutron star - main sequence star 12 pair favours stability against 
mass loss during the chaotic phase, given that results which apply to circular 
binaries can in some sense be carried over to highly eccentric binaries. One might 
in fact expect the stability criteria to be more severe, and this may explain the 
apparent dearth of cataclysmic variables in globular clusters (Shara 1996, but 
see also Grindlay 1996). It was recently suggested (Kochanek 1992) that most 
neutron star - main sequence star binaries formed by tidal capture in globular 
clusters would not survive the large angular momentum transfer involved in the 
circularization process. But this work assumed large tides (and hence large tidal 
angular momenta) all the way down to complete circularization. Given that many 
such binaries will spend the majority of their circularization time with rather small 
tides, tidal capture again becomes a viable mechanism for creating LMXBs. 

Some capture binaries will suffer a relatively minor chaotic phase. For ex- 
ample, given the rather low velocity dispersions found in open clusters and the 
relatively large cross section presented by a giant star, it is possible for such stars 
to capture other stars at "wide" initial separations. This may save such binaries 
from destruction and they may well go on to become interesting objects. On the 
other hand, giant captures in globular clusters are likely to be violent because in 
order for a capture to occur, the stars must approach each other extremely closely 
initially. 

It remains to be seen whether or not tidal capture binaries are important 
dynamically in globular clusters. Since many primordial binaries are destroyed 
in the core over the lifetime of a cluster, it may well be that, at least for very 
dense clusters, tidal capture binaries have an important role to play. The cluster 
calculations discussed above (Aarseth 1996) are presently being run on the HARP- 
2 machine at Cambridge with TV = 10,000, appropriate to large open clusters. 
Larger simulations with N = 25, 000 representing small globular clusters will soon 
be implemented on the HARP-3 machine which has recently been installed at 
Cambridge. 

Whatever the survival fraction of tidal capture binaries is, it is clear that the 
two-phase behaviour revealed by the chaos model for tidal interaction changes the 
way we now look at the process of tidal capture. 



12 The turnoff mass of a main sequence star in a globular cluster is around 0.8Mq. 
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